##############################################################
#Replication file for:
#The Importance of a Liberal Power's Attention to Democratic Elections Around the World
#Johannes Bubeck Ashrakat Elshehawy Nikolay Marinov Federico Nanni
##############################################################


#load data
library(readr)
library(haven)
rm(list = ls())

#read in data produced throygh 10_FifgureD.9_appendix.do
#change working directory
biasest <- read_dta("vdempolyarchy.dta")
##############################################################
##############################################################



#ggplot historgram
library(ggplot2)
plt0=ggplot(biasest, aes(x=bias)) + geom_histogram(binwidth=0.1) +
  scale_x_continuous(breaks=c(0,0.2,0.40000001, 0.60000002, 0.80000001,1),
                     name="bias") +theme_classic()+
  theme(axis.text=element_text(size=12),axis.title = element_text(size=13))
plt0 

##############################################################
##############################################################

#install needed libraries
#install.packages("Rmisc")
#install.packages("dplyr")
library(Rmisc)
library(reshape)
library(dplyr)


#our data - transformation needed through grouping using dplyr
dt<-biasest%>%
  dplyr::group_by(bias) %>%
  dplyr::summarize(mean=mean(t_vdems,na.rm=T),
                   lci=t.test(t_vdems,confl.level=0.95)$conf.int[1],
                   uci =t.test(t_vdems,confl.level=0.95)$conf.int[2])

#plot part 2 of the figure
pl1=ggplot(data = dt)
pl1 <- pl1 + geom_point(aes(x=bias, y=mean), color= "red", size = 3)
pl1=pl1+geom_errorbar(aes(x=bias,ymin=lci,ymax=uci,width=0.1)) +
  theme_classic()+
  scale_x_continuous(breaks=c(0,0.2,0.40000001, 0.60000002, 0.80000001,1),
                     name="bias")
#add text
pl1
pl1 <- pl1 + geom_text(aes(x=bias, y=lci,
                           label = round(lci,2)), size= 3.9, vjust = 1.2)
pl1 <- pl1 + geom_text(aes(x=bias, y=uci,
                           label = round(uci,2)), size= 3.9, vjust = -0.3)

pl1
pl1 <-pl1+theme(axis.title.x=element_blank(),
                axis.text.x=element_blank(),
                axis.ticks.x=element_blank(),
                axis.ticks.y=element_blank(),
                axis.text.y=element_blank(),
                axis.line=element_blank())
pl1=pl1 +scale_y_continuous(name="V-Dem Polyarchy") +
  theme(axis.text=element_text(size=12),axis.title = element_text(size=13))


pl1
##############################################################
##############################################################
#combine plots
#arrange
library(ggpubr)
theme_set(theme_pubr())
figure <- ggarrange(pl1,plt0,
                    ncol = 1, nrow = 2)
figure
setwd("/Figures")
ggsave("figured9.pdf")
##############################################################
##############################################################



